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We investigate the large N limit of pure M = 2 supersymmetric gauge theory with 
gauge group SU(A^) by using the exact low energy effective action. Typical one- 
complex dimensional sections of the moduli space parametrized by a global complex 
mass scale v display three qualitatively different regions depending on the ratio be- 
tween \v\ and the dynamically generated scale A. At large |f|/A, instantons are 
exponentially suppressed as A^ — oo. When \v\ ~ A, singularities due to massless 
dyons occur. They are densely distributed in rings of calculable thicknesses in the 
f-plane. At small |f |/A, instantons disintegrate into fractional instantons of charge 
1/{2N). These fractional instantons give non-trivial contributions to all orders of 
1/A^, unlike a planar diagrams expansion which generates a series in l/N"^, implying 
the presence of open strings. We have explicitly calculated the fractional instan- 
tons series in two representative examples, including the 1/A^ and 1/A^^ corrections. 
Our most interesting finding is that the 1/A^ expansion breaks down at singularities 
on the moduli space due to severe infrared divergencies, a fact that has remarkable 
consequences. 
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1 Introduction 



Quantum SU(A^) gauge theory has a single free parameter, the number of colours 
A^, or equivalently the number of interacting gluons A^^ — 1. Our best chance to un- 
derstand the strongly coupled dynamics of gluons is probably to perform an analysis 
of the theory at large |1|]. In the present work we will use exact results available 
for theories with eight supercharges (A/" = 2 supersymmetry) [0, ^ in order to try to 
better understand the large limit of four dimensional gauge theories. Some basic 
qualitative features were first discovered by the author while working on simple two 
dimensional toy models, which can be supersymmetric 0], but also non supersym- 
metric In four dimensions, we can unfortunately provide explicit calculations 

for M = 2 theories only. 

The paper is organized as follows. In Section 2, we give a general qualitative 
discussion of some of our main results, without entering into technical details. In 
Sections 3 and 4, we discuss in turn the structure of the moduli space at large A^, 
the nature and properties of the large A^ expansion of the low energy observables, 
and two full large A^ calculations of representative Seiberg-Witten period integrals. 
In Section 5 we discuss some open problems. 

2 General discussion 

Large N in SU (N) gauge theories and planar diagrams 

From the point of view of perturbation theory, the large A^ expansion is a reordering 
of Feynman diagrams with respect to their topology, with diagrams of genus h being 
proportional to l/N"^^. The original perturbative series has zero radius of convergence, 
but the large A^ reordering can produce convergent series in the gauge coupling con- 
stant at each order in l/N"^, thus providing a non-perturbative treatment. The series 
in l/N"^ is reminiscent of the genus expansion of an oriented closed string theory ^ 
of coupling constant k ~ 1/A^, and one may hope that a full non-perturbative equiv- 
alence between string theories and gauge theories could exist. This is supported by 
the empirical facts that hadrons are found on a Regge trajectory, and that quarks 
seem to be confined by the coUimation of Faraday flux lines in string-like structures. 
In a supersymmetric context, one can make this idea more precise by using the fact 
that D-branes, on which gauge theories live, can be viewed as solitons in closed string 
theories [^. 
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Large N and instantons 

A general amplitude in gauge theory can pick up contributions of different types. In 
addition to the terms having a Feynman diagram representation that are discussed 
in the context of the 't Hooft expansion, non-trivial field configurations invisible in 
Feynman diagrams can in principle contribute to the path integral. The most pop- 
ular configurations of this type are instantons |^ |10|. Instantons are responsible for 
important semi-classical effects like tunneling. In the case of A/" = 2 supersymmetric 
gauge theories, on which we will focus in the present paper, a non-renormalization 
theorem [|T^ implies that the low energy effective action S^s up to two derivatives 
terms picks up only a trivial one-loop term from Feynman diagrams, but has an in- 
finite series of instanton contributions. From the point of view of the large N limit, 
however, instantons are exponentially suppressed, and thus do not seem to play an 
interesting role [|12|. The instanton action is indeed proportional to in the 't Hooft 
scaling gyj^ oc 1/N. For Af = 2 super Yang-Mills, this seems to imply that Scs at 
large is simply given by the one-loop formula. One would have to look at higher 
derivative terms in order to obtain a non-trivial large expansion. 



More precisely, and as was emphasized in |jT2|, the effects of instantons of size 1/v 



are proportional in the one-loop approximation to 



e /^^YM = I -J ' (1^ 



where ^^ym is the gauge coupling constant at scale v, A the dynamically generated scale 
of the theory, and /3o oc A^ is given by the one-loop (3 function. The one-loop formula 
(ffl) is exact for TV = 2 super Yang-Mills, with Pq = 2N. It suggests that the only 



smooth limit of instanton contributions when A^ — > cxd is zero [|T2|. Large instantons 
(small v), if relevant, would produce catastrophic exponentially large contributions, 
and if one is willing to assume that the large A^ limit makes sense the only physically 
sensible conclusion is that instantons are irrelevant variables. In real-world QCD, 
instantons of all sizes can potentially contribute, and this led Witten to argue that 
the instanton gas must vanish |jl2|. In Higgs theories, like Af = 2 super Yang- Mills, 



the Higgs vevs introduce a natural cutoff v on the size of instantons, and for v large 
enough ("weak coupling") the instanton gas can exist but is just negligible at large 
A^. At small v ("strong coupling"), we run into the same difficulties as in QCD, and 
instantons must somehow disappear. 

Large N and fractional instantons 

An interesting puzzle is to try to understand the nature of the physics at strong 
coupling in the large A^ limit. Unlike in non-supersymmetric theories where one may 
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argue that the physics is dominated by planar diagrams, in A/" = 2 gauge theories 
we have already pointed out that some important observables do not have any non- 
trivial contribution from Feynman diagrams. In their derivation of the low energy 
effective action Ses, Seiberg and Witten made use of the electric/magnetic du- 
ality transformations of the low energy abelian theory. These transformations are 
very useful to understand the physics near singularities where magnetically charged 
particles become massless, but they say nothing about the large N limit at strong 
coupling. The purpose of the present paper is to elucidate this problem. We will 
see that instantons actually do not literally vanish, but disintegrate into fractional 
instantons of topological charge 1/{2N) (or a multiple of 1/(2A^)). Such fractional 
instantons give contributions of order 

-] =r' (2) 



V J 

and can obviously survive at large A^. The fact that fractional instantons could play a 
role in gauge theories has been suspected for some time, independently of the large A^ 
approximation. In the case of A/" = 1 supersymmetric gauge theories, chiral symmetry 
is broken by a gluino condensate (AA) = A^. Since Pq = 3N in this context, we see 
that the gluino condensate can be interpreted as coming from a fractional instanton 
of charge 1/A^ (a recent discussion of the gluino condensate can be found in [p!3| 



see 



also IjT^ for a discussion of the closely related problem of A/" = 1 superpotentials in 
confining vacua). In the following we will see that fractional instantons do play a 
fundamental role in A/" = 2 gauge theories as well, when the theory is analysed at 
large A^. 

Fractional instantons is a new type of contribution in the large A^ limit. Though 
it is known that Feynman diagrams generate a series in l/N"^ and that instantons 
must be exponentially suppressed, it is not obvious what is the analogous statement 
for fractional instantons. We will explicitly demonstrate below that in AT = 2 gauge 
theories, fractional instantons generate a series in 1/A^. The leading contribution can 
be of the same order of magnitude as the one coming from planar graphs. The first 
correction is then of order 1/A^, and would dominate any subleading correc- 
tions from Feynman diagrams. This proves directly from the field theory point of 
view that the string theory dual must contain open strings in addition to the familiar 
closed strings. This fact is actually consistent with our present knowledge on string 
duals to A/" = 2 gauge theories. In the supergravity approximation, the closed string 
background dual to such theories usually has unphysical naked singularities called 
"repulsons" [jl5i . String theory does actually make sense on such backgrounds thanks 



to the "enhangon" mechanism |16]: the constituent branes expand and form a finite 
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shell excising the singularity. The open strings responsible for the fractional instan- 
tons contributions must be open strings attached to these branes. Such open strings, 
and the fact that they would produce unusual corrections on the field theory 
side, do not seem to have been discussed in the literature. It would be extremely 
interesting if they could be used to match with non-trivial results of the type we are 
going to derive on the gauge theory side. 

Large N and singularities 

Possibly the most interesting result derived in this paper concerns the behaviour 
of the large N expansion at singularities on the moduli space. These results were 
anticipated in our study of two dimensional models 0, ||, |[. 

On the moduli space of vacua, the gauge group SU(A^) is broken down to U(l)^~^, 
and generically the low energy theory is a simple pure M = 2 abelian gauge theory. 
For some particular, "critical," values of the moduli, however, some additional hy- 
permultiplets can become massless 0. In the pure SU(A^) case on which we focus, 
such critical points can only be observed at strong coupling (the W bosons have then 
a mass of order A) and correspond to massless magnetically charged particles called 
dyons. Since those couple non-locally to the original photons, it is a priori not obvi- 
ous what the low energy theory can be, but we know from that the dyons couple 
locally to dual photons and generally produce a trivial free abelian theory in the IR. 
When the light particles are not mutually local with respect to each other, which 
can be achieved by adjusting several moduli to special values (higher critical point), 
a genuinely non-trivial CFT develops in the IR fl^. In gauge theories with quark 
hypermultiplets, very general types of critical behaviour can be obtained, see e. g. 

Commonly, trying to find a good approximation scheme to describe a non-trivial 
critical point can be subtle in field theory. A typical example is 0^ theory in dimension 
D. The theory has two parameters, the bare mass m (or "temperature") and the 
bare coupling constant g. By adjusting the temperature, we can go to a point where 
we have massless degrees of freedom, and a non-trivial Ising CFT in the IR. The 
difficulty is that the renormalized fixed point coupling g^, is large, and thus ordinary 
perturbation theory in g fails. It is meaningless to try to calculate universal quantities 
like critical exponents as power series in g, since those are (/-independent. Either 
the tree-level, ^^-independent contributions are exact and the corrections vanish (this 
occurs above the critical dimension, which is -De = 4 for the Ising model, and we have 
a trivial fixed point = well described by mean field theory), or the expansion 
parameter corresponds to a relevant operator and corrections to mean field theory 
are plagued by untamable IR divergencies. The way out of this problem is to use 
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expansion parameters on which the low energy CFT depends (corresponding in some 
sense to marginal or nearly marginal operators): we can use an e expansion by going 
to 4 — e dimensions, or a large N expansion by considering O(A^) vector models and 
O(A^) invariant critical points. 

The critical points on the moduli space of A/" = 2 super Yang-Mills are very similar 
to the Ising critical point below the critical dimension, and is very similar to the 
0^ coupling g. On any generic finite dimensional submanifold of the moduli space, 
one finds critical points, characterized by a set of critical exponents |19], that 



are independent of the number of colours of the gauge theory in which the CFT 
is embedded. These critical exponents cannot consistently be calculated in a 
expansion. Even the simple monopole critical points that are known to be trivial are 
not described consistently by the iV — oo hmit of the original gauge theory, because 
electric-magnetic duality is not implemented naturally in this approximation scheme. 
We will indeed explicitly demonstrate in Section 4 that, though the leading large 

approximation for the critical points is smooth ("mean field theory"), the 
corrections are IR divergent. 

Note that by adjusting a large number ~ of moduli, one can also study A^- 
dependent critical points. For these, the large A^ expansion can certainly be consis- 
tent. The A^ special vacua where a maximum number A^ — 1 of mutually local dyons 



become massless and which where studied in because of their relationship with 
Af = 1 vacua are of this type. 



3 The structure of moduli space at large N 
3.1 Brief review of = 2 super Yang-Mills 

The fields of pure SU(A^), Af = 2 supersymmetric gauge theory all transform in the 
adjoint representation of the gauge group and make up an TV = 2 vector multiplet 
which contains the gluons, gluinos and a complex scalar field $ with scalar potential 

l^=ltr[$,$t]2. (3) 
9 

The D-fiatness conditions [$, $"''] = are solved by 

($) = diag(ai, . . . ,aAr) (4) 

with 

TV 

^a, = 0. (5) 

i=l 
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Classically, the gauge group is generically broken to U(l)^^\ with W bosons of 
masses rriij = v^lflj — aj\. When some of the OjS coincide, the gauge symmetry is 
enhanced. Quantum mechanically, the moduli space is parametrized by the vevs of 
the gauge invariant operators 

Ui = {tT¥), 2<i<N, (6) 
or equivalently by the 0iS defined implicitly by the relations 

N N 

Uk = Y.^1, Y.^^ = ^. (7) 

i=l i=l 

The OjS are non-trivial functions of the 0jS, with ~ 0, only at weak coupling. 
In particular, it turns out to be impossible to have = a^, which implies that a 
non-abelian gauge group is never restored in the quantum theory. 

The low energy effective action is then generically a pure abelian U(l)^^^ theory. 
The leading terms in a derivative expansion can be written in A/" = 1 superspace as 



5'eff = — Im / (f'x 



where the {Ai, Wi)s are the massless abelian M = 2 vector multiplets satisfying 
= and T{A) is a holomorphic function called the prepotential. 

According to Gauss's law, the central charge Z of the supersymmetry algebra, 
and the masses 

Mbps = V2|^| (9) 

of the BPS states, only depend on the values of the fields at large distances, and thus 
they can be expressed in terms of JF only. Introducing the dual variables 

a-Di = (10) 
dai 

we have [§] 

N 

Z = ^{qiai + hia^i) (11) 

i=l 

where the integers g,, = 0; hi, J^i^i = 0; ^i^^ the electric and magnetic 

quantum numbers respectively. Low energy electric/magnetic duality is manifest on 
the formulas (|^) and ( pA]) which are Sp(2(A^ — 1), Z) invariant, and it turns out that 
the variables a and a£, are most naturally interpreted as sections of an Sp(2(A^ — 1), Z) 
bundle. 
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If one introduces the genus — 1 hyperelliptic curve C defined by 



N 



C : Y^ = J](X - (f),f - A^^ = P{Xy - A^^ (12) 



i=l 



equiped with a canonical basis of homology cycles (aj,/5j), and the differential form 

- 2mY' ^^^^ 

then a and a^) are simply given by 0, |^ 



ft: 



AT 



Asw, CLDi = f Asw- (14) 
It is natural to rewrite the equation for C as 

= (P(X) - A^)(P(X) + A^) = P+(X)P_(X) = J](X - X,+)(X - X,_) (15) 



i=l 



and to view the curve C as two copies of the complex plane with cuts joining the 
branching points Xj + and Xj _. The "electric" contour /3j defining encircles the 
cut from Xj + to Xj _. At weak coupling (small A), one then recovers Oj ~ 0j. For our 
purposes, however, it is not particularly useful to try to make a distinction between 
electric and magnetic variables. At strong coupling, the two notions are mixed by 
the non-trivial monodromies 0. The important point is that there is a one-to-one 
correspondence between possible BPS states and the integer homology of the curve 
C, in such a way that to any cycle 7 G Hi{C, Z) we can associate the central charge 

^(7)=/ Asw (16) 

Singularities on the moduli space are obtained when Hi{C,Z) degenerates, or equiv- 
alently when two of the branching points coincide. Since the polynomials P+ and P_ 
obviously have no common roots, this can happen only when the discriminants of P+ 
or P_ independently vanish, 

A(P+) = n(^.,+ -^.,+)' = or A(P_) = n(^^,--^.,-)' = 0. (17) 

i<j i<j 
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3.2 The large N limit 



3.2.1 Generalities 

At large A^, it is convenient to parametrize the N — 1 dimensional moduli space by 
the density function 

1 ^ 

i=l 

which satisfies 

%,v(0)rfV = l, /#^(0)rfV = O. (19) 



At weak coupling, Pn{<P) simply gives the density of complex eigenvalues of the Higgs 
field. The differences 0j — (pj are related to the masses of the W bosons. If the lightest 
W has a mass of order mw, then the heaviest W will generically have a mass of order 
Mw ~ y/NrriY/. Particles having a mass growing with A^ are problematic in the large 
A^ expansion. This is due to the fact that the amplitude for the propagation of a state 
of mass M for a time t involves the factor exp{—iMt) (this argument was used long 
ago by Witten for baryons |21|). A consistent large A^ limit, for which propagators 
of W bosons have a smooth expansion, can thus be achieved only if the heaviest Ws 
have a mass of order N^, and thus the lightest Ws have generically a mass of order 
1 / \/N. For one dimensional distributions corresponding to the case where all the 0jS 
are aligned, the lightest Ws have a mass of order 1/A^. This correct large A^ scaling 
was already considered in [|^. Under the above conditions, the distribution Pn{<P) 



will typically have an A^ — oo limit p(0) which is the sum of a smooth function with 
compact support plus S function terms. We will consider examples below. 

A good strategy to study the large N, N — 1 dimensional, moduli space, is to 
consider well-chosen finite dimensional subspaces. The simplest subspaces Aip are 
one-complex dimensional and are parametrized by a complex mass scale v. One picks 
up a particular fixed distribution 0j characterized by p(0) and considers = vipi 
or equivalently 

p,(0) = lp(^)- (20) 

When p(0) is a smooth function, then the limit — oo is a weak coupling limit 
since the lightest Ws have a mass proportional to v. If p(0) has some 6 function 
singularities, which means that a large number (of order A^) of the 0jS coincide, then 
we are at strong coupling for any v. The study of the simple subspaces M.p can 
reveal a great deal of physics, including the transition from weak coupling to strong 
coupling and the presence of the simplest singularities. To obtain more complicated 
singularities {kih. order critical points), one has to consider k dimensional subspaces. 
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or introduce quark flavors and quark mass parameters. For the purposes of tlie present 
paper, it is enougfi to stick to tlie one- dimensional subspaces M.p defined by (|^). It 
is then natural to introduce the dimensionless ratio 



r = I (21) 



and the polynomials 



N 



pix) = Hix - 0,) , p±(x) = pix) T (22) 

i=l 

and to rewrite the fundamental equations (|12D, ([T3|) and (0) in terms of dimensionless 
variables, 

C: =p{xf =p+{x)p4x), (23) 

X = Xsw / v = xdp / {2iTTy), (24) 

z{-f) = Z{^) / V = (fx. (25) 
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3.2.2 Singularities 

The locus of singularities S corresponds to the points where the discriminants (|T^ 
vanish, or equivalently when one of the polynomials p+ or p_ has a multiple root. 
The discriminants are themselves polynomials of degree — 1 in r^, and thus we 
will have generically 2N{N — 1) singularities where a dyon becomes massless on Aip. 
These singularities are arranged on 2{N — 1) concentric circles. To better understand 
the structure of iS, it is useful to study first the large N distribution of the roots of 
the polynomials p±. The equation p±{x) = can be rewritten 

In \ r^p(x) I + i arg (±r^p(x)) = 0. (26) 

The real part of the equation is of order while the imaginary part is of order iV°. 
In the leading large approximation, both p+ = and p_ = thus reduce to 

E;v(x) =^b(x)|i/^ = l/|r|. (27) 
The N oo limit of the "envelope" En is easily evaluated, 

E^{x) =exp Jd^(l)p{(j))\n\(f)-x\. (28) 

When X is in the support of p, the above approximate formula only gives an average 
value of Ei\f, which in fact has very sharp wells around the 0jS, Ei\f{(j)i) = 0. In the 
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Figure 1: Plot of the function En=25{x) corresponding to the distribution ( ^2]) with 
(pi = 2{i — 1)/N — 1 {1 < i < N) (thin hne) together with the envelope E^o (thick 
line) . 



large expansion, the wells are approximated by infinitely thin lines joining the 
smooth envelope (see Figure 1 for an example). Any correction to this picture 
would be exponentially small at large N, an instanton effect. At large 



E^,(x) 



\x\ 



[l + 0{l/\x\)) 



(29) 



Equation ( jn\ ) is now easily solved. When l/|r| < min^-gc -E'oo(a;), the plane z = l/|r| 
can only intersect z = En along the thin wells: the roots of p+ and p_ coincide with 
the classical roots 0^. This corresponds to the instanton dominated region, and those 
are suppressed at large A^. When \r\ decreases and l/|r| > min^jgc -£"00(2^), z = l/|r| 
and z = En continues to meet along some of the thin wells, but the roots nearest 
to the minimum of E^o are now arranged along an inflating curve. This curve gives 
the general shape of the "enhangon" discussed in [0 . Note that the enhangon exists 
for all r if the distribution p has some 5 function singularities. It has in general 
several connected components associated with the different connected components of 
the support of p. At large all the roots are located on the enhangon whose 

various connected components eventually merge to form asymptotically a circle of 
radius l/|r| according to (p^). 



For example, let us consider the simple uniform distribution 



Po(0) 



l/vr if 
if \<h\ 



01 <1 
> 1. 



(30) 
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It yields an envelope 



For |r| > i/e, all the roots are classical. At |r| = ^/e, the enhangon begins to inflate. 
It is a circle of radius a/1 — In |rp. At |r| = 1, the enhangon has eaten up all the 
classical roots, and for |r| < 1 it is simply a circle of radius l/|r|. For the purpose 
of illustration, we have depicted in Figure 1 the case of a unidimensional uniform 
distribution 

^^^'^^"loif |0|>1 

for which 

r (1 - x)(i-")/2(_i _ x)(i+")/Ve if X < -1 
EocAx G M) = < (1 + x)(i+^)/2(i _ 2,)(i-x)/2/g - 1< X < 1. (33) 

[ (l + x)(l+")/2(_i +2,)(l-x)/2/g if^>l. 

Now that we understand the location of the roots of p±, it is straightforward 
to deduce the singularity locus. Generically, roots are either equal to the classical 
values, or are smoothly distributed along the enhangon. Singularities can then occur 
only when the classical roots are eaten up by the inflating enhangon, or when two 
connected components of the enhangon merge. This is particularly transparent on 
Figure 1, where real roots join the enhangon by pairs and then become complex. 
We thus deduce that the singularities are densely distributed on rings in the r-plane 
corresponding to the melting of classical roots with the enhangon (singularities of the 
first class) or with the melting of several connected components of the enhangon with 
each other (singularities of the second class). In Figure 2 we have represented the 
space Aipo with the associated singularity locus S (only singularities of the first class 
are present in this example). We will see in Section 4 that the large expansion 
is afflicted with power-like divergences near singularities of the second class, while 
logarithmic divergencies also mix up near singularities of the first class. 



3.2.3 Special cases 



We have discussed above how the singularity locus S can be deduced for a given 
density p. It is amusing to study the inverse question: can we find a density p 
that would correspond to an a priori given locus S7 This amounts to inversing the 
functional relation ( pSj) between E^o and p. If we limit ourselves to the consideration of 
unidimensional distributions (by constraining all the 0jS to be real), the mathematical 
problem is very similar to the one encountered in the study of the leading N ^ oo 
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r-plane 




^1 singularities 

Figure 2: The one-complex dimensional subspace Aip^ corresponding to the uniform 
distribution ([30| ) with the associated ring of singularity S. 



limit of simple matrix integrals p2[. If we introduce 



uj{x)= ld<f)-^, (34) 
X — d) 



we have 



and dl) yields 



p(0) = ^ (^(0 + - ^(0 - ^e)) , (35) 
uj{(t) + ie) + - ie)] . (36) 



2ti 

J_dE^^ 1 

As an application, let us try to find a simple unidimensional smooth density p with 
a connected compact support [—1, 1] for which all the singularities are distributed on 
a single circle (a singularity ring of zero thickness). From the discussion of Section 
3.2.2 it is clear that this can only happen when dE^o/dx = on the support of p. 
The only solution to ( pG]) consistent with uj{x) ~ 1/x at infinity is then 

U!(x) = ; (37) 

v/(^-i)(x+i) 

and ( p5|) implies 

P(0) = ^i=^- (38) 



This is the density distribution of the roots of the Chebyshev polynomials T^, for 
which (pi = cos(7r(l/2 + i)/N). This very special distribution was studied in [pO| . 
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Equation (^) then gives the radius of the singularity circle |r| = 2. Clearly, many 
generalisations could be studied. 



3.2.4 Corrections to iV = oo 



The singularity rings of the first class separate an outer purely instanton-dominated 
region from an inner strongly coupled region which cannot be described by semi- 
classical physics (the outer region only exists when p does not have any 5 function 
singularities). In the outer region, corrections at large N are exponentially suppressed. 
To the contrary, we can a priori expect non-trivial corrections in the inner, strongly 
coupled, regions of moduli space. We will elucidate the nature and the main properties 
of these corrections in Section 4, but we can have a first instructive look at them by 
investigating the envelope or equivalently the shape of the enhangon on a simple 
example. Let us choose (pi = 2{i — 1)/N — 1 as in Figure 1, and let us calculate the 
corrections to (|33D for x > 1. By applying Euler's formula 

° g /(a + k'-^) = ff - - °>'^f - ^'°» + Odin-) (39) 



n 



k=0 



we obtain 



Ejv,i(x) = E^,,ix) (l - t4 In ^ + 0(l/iV2)) for a; > 1. (40) 

This formula displays two of the main features of the corrections to iV = oo: they are 
of order {not of order l/N"^) and they diverge near the singularity locus which 



corresponds here to |r| (l/2)e^ or equivalently a; ^ 1^ (see (^) and (|33D ). A 
third feature, that the corrections are given by series of fractional instantons, is also 
to be expected if one consider that (^) depends on 1/r. 



4 Calculation of the and corrections 

This Section is devoted to a detailed study of the corrections to the leading N ^ oo 
approximation in two particular, but generic, examples, encompassing the cases of 
both types of singularities (first class and second class) . Since we want to consider the 
purely strongly coupled region, where the physics cannot be described by instantons, 
we can put most of the singularities at infinity. Yet, we will keep one singular circle 
at finite distance, so that we can discuss the physics associated with singularities. 
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-1 -0.5 0.5 1 1.5 

Figure 3: Plot of the function E]\r=25{x) corresponding to the distribution (|ll| ) (thin 
hne) together with the envelope E^o (thick line). The roots x_ and Xi of the poly- 
nomial p_ discussed in the text (equations ( ^^ and (^S])) are also indicated for a 
particular value of 1/r. 

4.1 An example with a first class singularity 

We choose the distribution pat to be 

Piv(0) = 5(0 + l/N) + ^ 5(0 + l/N - 1). (41) 

The associated function Ej^ and envelope E^o are depicted in Figure 3. The parameter 
r will be real and positive, without loss of generality (the final formulas can be trivially 
analytically continued). We have a singularity of the first class corresponding to the 
melting of the classical root Xd = 1 — l/N with the enhangon. 

4.1.1 Infrared divergences and fractional instantons 

A very simple, yet very instructive, calculation one can do is the critical value of 
r for which the singularity occurs. The leading ^ oo approximation = 1 can 
be read off immediadely from Figure 3. The exact value can be found by noting that 
it corresponds to the point where the two positive real roots of the polynomial p_ 
defined by ( ^2]) coincide, 

_ N 

~ (iv- l)i-W' ^ 
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The large N expansion of the exact formula ( ^21) has strange looking logarithmic 
terms, 

In + 1 

rc = l + ^^ + o(l/iV). (43) 

A standard, well-behaved, large N expansion cannot produce such (lnA^)/A^ correc- 
tions. Remarkably, the same kind of (lnA^)/A^ corrections to the critical parameter 
were found recently by the author in a study of a two dimensional QFT that was 
argued to be very similar to gauge theories with Higgs fields |^. These (In N)/N 
corrections were the consequence of the breakdown of the large N expansion near the 
critical point due to infrared divergences (a full calculation can be found in Appendix 
C of 01). We will see below that the qualitative physics in four dimensions is strictly 
similar to the physics previously discussed in two dimensions @, |^, ^. 

The next simple calculation one can obviously do is the large N expansion of the 
roots of the polynomials p+ and p_. These roots also enter explicitly in the Seiberg- 
Witten period integrals (|25|). We will focus on the two positive real roots of the 
polynomial p- which exists for r > and coincide at r = rc. One of these two roots 
that we call x_ is simply 

x. = l- l/N -l/r^ + 0{l/r^^). (44) 

The non-trivial corrections to the classical root Xc\ = 1 — 1/N are exponentially small. 
The other root Xi has a non-trivial and very interesting large expansion that can 
be straightforwardly deduced from (p6|). 



1 r -|- ln(r — 1) 1 , ,,9 2rln(r — 1)\ , , , 

^1 = + ^ - 1))' - + 

This very simple formula displays all the features that were advertised in Section 2. 
In particular, by writing ln(r — 1) = Inr + ln(l — 1/r) = Inr — X]a;>i ^/{kr'^), we 
see that each order in 1/N is given by a series of fractional instantons of fractional 
charge 1/{2N), mixed with the one-loop diagram contribution in Inr. At r = 1, 
the corrections are blowing up like the logarithm of the mass of the lightest degrees 
of freedom, an infrared divergence. Using (|^) and (|^) to solve x_ = Xi, one can 
recover the expansion (|43D with the (lnA^)/A^ correction. 

Of course, roots like Xi are not directly observables. The physics is encoded in 
the Seiberg-Witten periods (pSf), and though they depend on the roots, one might 
argue that cancellations could occur and that the final result could be smooth. It is 
also conceivable that only corrections in l/N"^ could show up in the final results. To 
answer rigorously these questions, we will compute in the next subsection the Seiberg- 
Witten period integral z corresponding to the cycle 7 which vanishes at r = Tc. It 
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turns out that the only potential divergences come from P5| ) in this case, and that 
the quahtative features of the expansion for xi and for z are the same. 



4.1.2 The calculation of the period 

We thus proceed to compute the large expansion of 



^ = A = - r'dx , "^^'^"^^ (46) 



7 



with 

p{x) = {x + l/N)^-\x-l + l/N). (47) 
It is convenient to trade the variable x for u = l/{r{x + 1/A^)), in terms of which 

p du (1 - ru/N){l - (1 - 1/N)ru) 



mr 7„_ u2 v^(l - ™)2 - -aiv 



(4J 



with 



u_ = ^ (^l + l/r^ + C(l/r2^)), (49) 
«i = 1 + ^^^^^ + ^ ((Mr - 1))' + ^ ln(r - 1)) + 0(l/iV^ (50) 

At large A^, it is tempting to neglect the term in (^8]), which is exponentially 
small on most of the integration region. This is correct except for m ~ mi ^ 1, when 
u"^^ is not negligible, and for n ~ n_ 1/r, when {ru — 1)^ can be smaller than m^^. 
We will thus distinguish three different contributions to z, 

z = zo + z^ + Zy , (51) 

with 

,„.jLrf!^(i-i^)(i_(i^^) 1 (52) 

tnr M^V A'/V A* / ru — 1 

ZTrr^ AT A AT 7 ^(^^ _ i)2 _ ^2iv ru-l/^ ' 

i7crj^_u^\ nJ\ N / I ^(to - 1)2 _ w27V ru-lM ^ 
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where U e]u-,ui[ is an arbitrary constant, for example U = (m_ + mi)/2. The integral 
zq can be explicitly evaluated, and its expansion is given by 

Finding the asymptotic expansion of z-^ at large iV is a bit more subtle. The integral 
would be exponentially small if not for the integration region near u\. The term v?^ 
cannot be neglected in a region of size XjN around U\. The idea is to use a new 
variable w defined by 

u — u\ — w/N. (56) 

The large N expansion of the integrand can then be straightforwardly deduced, and 
one is left with several definite integrals on the variable w e [0, A^(wi — U)\. The 
integration region can be extended to [0, oo[ at the expense of neglecting exponentially 
suppressed terms. The calculation of the and l/N"^ corrections involve five non- 
trivial definite integrals that can be explicitly evaluated. Instead of writing down 
all these complicated formulas, we have illustrated the procedure on a simpler, but 
qualitatively similar, example in the Appendix. Fortunately, the final result is simple- 
looking, 

i-KZ^ In 2 1 ( (, , 1 1 1 o 



+ ( fln(r - 1) +r+ -ln2Vn2 
N Nr \\ ^ ' 2 J 24 



(57) 



The evaluation of z^ proceeds along the same lines, except that now the relevant 
integration region around u = is exponentially small. The asymptotic expansion 
to all orders in 1/N can then be obtained in this case by changing the variable to w 
defined hy u — + w/r^'^^. We get 

^^ = ^ln2 + 0(l/r''). (58) 

Putting all the contributions together, we finally end up with 

mz r — 1 , 1/, r — (r — 1) ln(r — 1) \ 

l^-—r + ^'''+N [-^-r+-^ln2+'- L-A Lj (59) 

^ (ln(r - 1))^ + (In 2) ln(r - 1) + ^ (ln2)2 - ^] + 0{1/N% 



N^r ^ // ' V /V ^ 2 ' ' 24 



This final formula displays all the remarkable properties of the large N expansion 
of jV = 2, SU(A'") super Yang-Mills theory in four dimensions that we have already 
discussed at lenght in the previous Sections. We note that genuine divergencies when 
r ^ 1 show up only at order but non- analytic terms are already present at 

order 
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r = 1 .02 



r= 1 
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Figure 4: The enhangon for the distribution ( |60D at r = 1.02, r = 1, r = 0.98 and 
r = 0.5. The dots correspond to the actual roots of the polynomial — 1/r^^ 

for N = 30. At r = = 1, two roots coincide, the two components of the enhangon 
merge, and we have a second class singularity. At small r the enhangon is a circle of 
radius 1/r. 

4.2 An example with a second class singularity 

Let us now pick to be even for convenience and consider the case of the distribution 

PN{<P) = l{5i<P-l) + S{<P+l)) (60) 



which corresponds to the polynomial p{x) 



[X 



1)^/2 (ggg The enhangon 



has two components at large |r| which merge at a second class singularity for |r| = 1. 
The shape of the enhangon, together with the roots of p{x)'^ — Xjr^^ are depicted in 
Figure 4. 

If we take r to be real for concreteness, then the critical value of r is exactly given 



by 

r-c = 1 

and corresponds to the merging of the branching points 



X-t = ±v^l — l/r^. 



(61) 
(62) 



The formulas (^T]) and (|62D are to be compared with (|4^j43|j4^j45|) . They exhibit 
one important difference between first class and second class singularities: the large 
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N expansion of the roots (or the shape of the enhangon) is perfectly smooth in 
the case of second class singularities. However, for the same general reasons as in 
the case of first class singularities, the large N expansion does break down as well. 
To demonstrate this, we have to consider the Seiberg-Witten periods which are the 
physical observables. 



1 r+ , xp'{x) N r du l-l/{ur^) 



m J,_ - l/r2^ «7r //,2 u \l 1 - ' ^^^^ 

where u = l/{r'^{l — x"^)). A straightforward calculation along the lines of the pre- 
ceding subsection or of the Appendix then yields the large expansion 



— = - l/r^ - 2 In [I + - l/r^j - 2 In r - -j^ 

We get non-analytic terms at leading order, and genuine divergences at order l/N"^. 



5 Open problems 

We have seen that, despite the fact that the low energy effective action oi M = 2 
super Yang-Mills is, in some sense, generated by instantons only, it has a surprisingly 
rich and interesting large A^ expansion. 

Low- dimensional sections of the moduli space can be easily analysed. For example. 
Figure 2 is very reminiscent of the structure of the moduli space for N = 2. This 
suggests that the spectrum of BPS states and the curves of marginal stability could 
have a simple description at large A^, and that the methods of |^ may be useful. 

We have elucidated the fate of large instantons at strong coupling and large A^ 
in A/" = 2 super Yang-Mills: they disintegrate into fractional instantons. These 
fractional instantons are responsible for a new class of contributions in the large A^ 
expansion, generating a series in 1/A^ in addition to the standard series in l/N"^ from 
Feynman diagrams. These new contributions were brought to the fore particularly 
clearly because we have studied observables which pick up a single one-loop term 
from Feynman diagrams. More general observables will probably have non-trivial 
contributions from both types of terms. Note that the leading corrections to A^ = oo 
are dominated by fractional instantons. It would be very interesting to find direct 
general arguments explaining why fractional instantons contribute a well-behaved 



20 



asymptotic series in The fate of instantons in real world QCD remains of course 
an open problem, but it is fascinating to contemplate the possibility that phenomena 
qualitatively similar to the one studied in the present paper could occur. This is 
suggested by the fact that the basic difficulty with instantons — IR problems due to 
the Landau pole for large instantons — are independent of supersymmetry. 

The standard 't Hooft analysis of perturbation theory at large N in SU(A^) gauge 
theories |T[] is consistent with a dual description in terms of closed oriented strings. 
Our analysis of the fractional instanton contributions shows that open strings must 
also be present. Though adding open strings to a closed string theory is a non-trivial 
step, it is satisfying that fractional instanton terms can be interpreted in a string 
picture at all. 

The breakdown of the large expansion at singularities may come as a surprise 
(it was anticipated in the recent studies of two-dimensional models by the author 
[0, 1^, H), but we have argued that the physics from the field theory point of view is 
the same as the one that produces IR divergencies in corrections to mean field theory 
below the critical dimension. On the other hand, the interpretation from the string 
theory point of view is much less clear. A naive application of the UV/IR relation 
seems to indicate that the string dual is not UV renormalizable, but this would be 
a rather strange property for a string theory. Another possibility, which was pointed 
out to me by I. Klebanov, is that the divergencies are due to tensionless strings at the 
singularities. This interpretation requires the appearance of divergences (or at least 
non-analytic contributions) at leading order, which is not the case in (|59D. This might 
nevertheless be made consistent by correctly identifying the terms generated by closed 
strings or open strings diagrams. This point clearly deserves further investigation. 

Another fascinating line of research, suggested by the divergencies of the 1/A^ 
expansion near singularities, is that it might be possible to extract finite universal 
string amplitudes by taking N —* oo and r — ^ rc in a correlated way, along the lines 
of [p^. A preliminary investigation will be presented in a forthcoming paper p5]. 
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Appendix: a simple toy integral 



Let us consider the integral 



By changing variable to y = , In is straightforwardly expressed in terms of Euler 
beta function, 

In = —B(l/N,l/2)-l = -1- (66) 

^ / ' / ^ ^ r(l/2 + l/A^) ^ ' 

The large N expansion can then be readily found, 

2 In 2 2(ln2)2-7rV6 , 
^iv = ^ + -^-^2 + 0{l/N'). (67) 

The difficulty one encounters to understand the large asymptotic expansion directly 
on the integral formula (^) is the same as for the integral (0) studied in the main 
text: the integrand is exponentially small (an instanton effect) except for a region 
of size near x = 1. One then use the variable y defined by x = 1 — y/N, and 
expand the integrand, which yields 

!/■",/ 1 X If", y^e-» 



The integration region can be extended to y E [0,oo[ by neglecting terms of order 
. The integral at order is elementary and gives the (21n2)/A^ correction. 
The integral at order l/N"^ is of the type of the integrals encountered in the main 
text. By going to the variable y' = l/vT~^-"e~^ and integrating by part, it can be 
related to the dilogarithm function Li2 and thus explicitly evaluated. 
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